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Abstract 

The spectral determinant of the Schrodinger operator (— A+V(x)) on a graph is computed 
for general boundary conditions. (A is the Laplacian and V(x) is some potential defined on 
the graph). Applications to restricted random walks on graphs are discussed. 



1 Introduction 

The study of spectral properties of the Laplacian operator on finite graphs began about fifty 
years ago. Many different domains are interested in the knowledge of those properties - let us 
simply mention organic molecules Q, superconducting networks vibrational properties of 
fractal structures ||], weakly disordered systems [§] and, more recently, quantum chaos ||. Of 
course, mathematicians § are also interested in that subject. 

Let us come back, for the moment, to the physics of disordered systems. In 0], the authors 
emphasized the central role played by the spectral determinant of the Laplacian in the compu- 
tation of the weak localization corrections. By constructing the Green's function on the graph, 
they obtained a compact form for this determinant. (See also || where a path integral approach 
is developed; in particular, a trace formula for the Laplacian on a graph || is recovered). 

Recently, the result of 0] was generalized || [h| to the spectral determinant det(H + 7) (= 
5(7)) with H = —A + V(x). V(x) is some external potential defined in each point x of the 



graph and 7 is a constant (spectral parameter). In |1C], the computation was done with the 



help of a path integral representation of the spectral determinant and also using time-dependent 



harmonic oscillator properties. Schrodinger operators have also been considered in 11] where the 
scattering matrix is computed for graphs made of one-dimensional wires connected to external 
leads. 

All this was done assuming continuity of the eigenfunctions at each vertex. 
Nevertheless, this "natural" (!) assumption is highly questionable. 

For instance, in 0|, the authors argue that the reduction of a realistic system of cou- 
pled tubes to a graph model is far from being obvious. In particular, serious problems arise 
from the finite thickness of the tubes, the geometry of the connection regions and also from 
eventually applied external fields. Analyzing in details a model of junction (what they call the 
"geometric-scatterer junction"), they suggest that it would be more appropriate to consider 
general boundary conditions on the resulting graph. 
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This is the point of view we will take up in this paper when computing the spectral determi- 
nant]]. Moreover, we will show that "playing" with the boundary conditions allows us to study 
some properties of closed random walks on any graph. For instance, it is possible to count the 
number of such walks when the number of backtrackings on each of them is fixed ]lq]-[o]. 



The paper is organized as follows. In section ||, we set up the notations that will be used 
throughout the paper. Section ||| will be devoted to the computation of the spectral determinant 
5(7) for general boundary conditions. Another expression for 5(7) will be derived in section ||] 
for permutation-invariant conditions. Applications to countings of restricted random walks on 
any graph will be discussed at the end of this section. Finally, a short conclusion will be given 
in section |5[ 

2 Definitions and notations 

We consider a graph Q made of V vertices, numbered from 1 to V, linked by B bonds of finite 
lengths. The coordination of vertex a is m a (X)o=i m a = 2B). 

On each bond [a/3], of length l a p, we define the coordinate x a p that runs from (vertex a) 
to l a /3 (vertex f3). We will also use xp a = l a p — x a p. 

To avoid cumbersome notations, $ being some function defined on the graph, we will simply 

write J[af3] $ for /o 1 " 5 dx af3 • 

An arc (a/3) is defined as the oriented bond from a to (3. Each bond [a/3] is therefore 
associated with two arcs (a/3) and (/3a). In the sequel, we will consider the following ordering 
of the 2B arcs: (lai)(la 2 ) . . . (la mi )(2/3i) . . . (2/3 m2 ) . . . 

Concerning the eigenfunctions <f> of H on the graph, we define on each link [a/3]: 

4>{al3) = I™ <t>{x a p ) ; 4>(/3a) = Hw. 4>(xf3 a ) (1) 

{^'{ap) ls ^ ne outgoing derivative at vertex a along the arc (a/3)). 

For the Green's function G(x,y) (x G [a/3], y anywhere on the graph), we similarly define: 

G(af3)(y)= Hm G(x a/3 ,y) ; Gm a) {y)= lim G(xp a ,y) (3) 
n i ( \ y dG(x Q/3 ,y) , dG(x/3a,y) 

With the above quantities, we can build the four (2B x 1) vectors </>, <//, G(y) and G'(y), 
respectively of components (q/3) , <// (a/3) , G (a ^{y) and G' {a(3) (y). 

In those conditions, the generalized boundary conditions for the operator H on the graph 
can be written: 

C + D <// = (5) 

1 Spectral properties of graphs with general - even random - boundary conditions imposed at the vertices have 
already been studied in the context of quantum chaos 
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where C and D are two (2B x 2B) constant matrices that don't depend on 7. 



In [18], the authors established the conditions for the operator H to be self-adjoint: CD + 
must be self-adjoint and the {2B x 4B) matrix (C, D) must have maximal rank 2B. 

Local boundary conditions connect, for each vertex a, the 4>(aPi) ,s to the <^( a/ g.)'s, i,j = 
1, . . . , m a . For such conditions, C and D can be chosen block-diagonal, the square block C a (or 
D a ) being of dimension m a . If, in addition, we assume that, for each vertex a, the conditions 
are invariant in any permutation of the nearest neighbours of a, we can write: 

C a = c a l + t a F a (6) 
D a = d a l + w a F a (7) 

where 1 is the unit matrix and F a is a matrix with all its elements equal to 1. The constants 
c a , d a , t a and w a characterize the boundary conditions in a. 

Remark that c a and d a can't both vanish because of the maximal rank condition^]. 

It is easy to realize that the quantity c a 0( a ^) + d a 4>^ aj3 ^ doesn't depend on i, i.e. it is the 
same for all the arcs starting at a. To conclude this section, let us mention the two limiting 
cases: 

i) d a = that ensures the continuity of 4> at vertex a (thus 4>{a) is defined) and leads to 

j=i ^ a ' 

(A a = corresponds to Neumann boundary conditions). 

ii) c a = 0. In that case, all the outgoing derivatives in a are equal. 

Now, we turn to the computation of the spectral determinant 5(7) (= det(H + 7)) of the 
operator H = —A + V(x) defined on the graph with boundary conditions given by (|5|). 

3 General boundary conditions 

As in 0, we construct the Green's function G(x,y) on the graph: 

(j + H)G(x,y) = 5(x-y) (8) 

and use the relationship: 

[ G(x,x)dx = d 7 lndet(# + 7) (9) 

J Graph 

In this section, we will consider, for each bond [af3], two independent solutions, ip a p and 
ipp a , of the equation: 

2 Moreover, this condition imposes that, at least, one of the two matrices, C a or D a , is invertible and can be 
set equal to 1 because of the homogeneity of condition (^). Finally, self-adjointness of CD + implies that only two 
real parameters are, actually, necessary to characterize the boundary conditions at each vertex a |isj ]. This will 



appear explicitly in section 1.2 where our results are expressed in terms of the two parameters r] a and p a defined 
in eqs.M H- 
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(H + 7) V = 
Those functions are chosen to satisfy: 

4> a p(a) = 1 ; ipapiP) = 



Their wronskian may be presented as: 

dV>/3a , Wap _ dipp c 



W a p = 1p a p 



dx 



a/3 



4>Pc 



dx a/3 dx a/3 



-(a) 



dip a p 



dx 



a/3 



(P) = Wp c 



We also define: 



"5 



dx^ 



(xp a = 0) 



So, let us construct this Green's function G(x,y). We first suppose that y belongs to 
link [ab]. 

If x is located on another bond [a/3], we have: 

G{x,y) = G {af3) (y)^ a p{x) + G^ a ){y) tpp a (x) 
Taking the derivative in a and using ( p^Jl^) , we get: 

G '(ap)(y) = G (aP){v) 1p'ap( a ) + G (Pa){y) W a p 

On the other hand, if x belongs to the same bond [ab] as y, G(x, y) must satisfy, when e 



G(y-e,y) = G(y + e,y) 



dG 

dx 



x=y-e 



dG 

dx 



+ 1 



x=y+e 



This leads to: 



x <y G(x, y) = G {ab) (y) ^ a b{x) + G {ba) (y) ip ba (x) + 

x > y G(x, y) = G (ab ) (y) ip a b(x) + G {ba ) (y) ip ba (x) + 

{x < y means that point x is closer to a than y). 
For the derivative in a, we obtain: 

G '(ab)(y) = G (ab)(y) ^ab(a) + G (ba)(y) W ab + ^ ab (y) 

Eqs. (16,^) can be written in matrix form: 

G / (y) = NG(y)-L(y) 
where N is a (2B x 2B) square matrix with elements: 

N( a p)( m ) = ^a^p^'apiot) + 5 ar) 5p^W a p 



W ab 
Tpba(y)^ab{x) 

w ab 



and L(y) is a (2B x 1) vector: 

L (y)(a/3) = ~ {SaafybTpabiy) + S ab Sp a 1p ba (y)) (24) 

Eq.(22) together with the boundary condition (fST) lead to: 



We deduce: 



G(y) = TL(y) (25) 
with the square matrix T = (C + DNy 1 D (26) 



G(ab){u) — T(ab)(ab)L(y)(ab) +T(ab)(ba)L(y)(ba) ( 27 ) 

G(ba){y) = T(ba)(ab) L (y)(ab) +T(ba)(ba)L(y)(ba) ( 28 ) 

and, after simple manipulations: 

G(y,y) =T (a6)(o6 ) (-^(y)) + T( 6a)(6a ) (-Vla(y)) + 

+ (- T (a6)(6a) " r (6a)(a6) + ^) (Mv) ^ba{y)) (29) 

To take the trace of G, we must first integrate ipab(y)i rfbaiv) an d Tpab(y)ifJba{y) on [aft]. We 
have shown in Q that: 

I y ah = -d^' ah (a) (30) 

/ ipl = -d^' ba (b) (3i) 

J[ab] 

lpab4>ba = -d-yWab (32) 



ah] 



Thus: 



G(y, y) = T^ ab )(ab)d^' ab (a) + T iba ^ ba) d^' ba (b) + (r (a6 )( 6o ) + T (ba)(afe ) ) SyWab - 

ab] v 7 

' (33) 



W ab 



Now, we sum over all the bonds. With the definitions of the matrices ./V and T, we obtain 
(C and D don't depend on 7): 



/ G(y,y)=Tv((C + DN)- 1 d 1 (C + DN))-d 1 In WaJ 
jGra P h \H / 



(34) 



Finally, with the observation that Tr((C + DN)- l d-y{C + DN)) = d 7 lndet(C + DN), we 
get the spectral determinant (up to an inessential multiplicative constant): 

5( 7 ) = det(H + 7) = TT TIT" det (C + DN ) ( 35 ) 

[a/3] W "P 

where II [a/3] • • • means product over all the bonds. 
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The expression (35) is valid for quite general (even non-local) boundary conditions)^. 



It is worthwhile to mention that eq. (|35|) can be, heuristically, recovered by a path integral 
approach (see (p!o|]) for the Neumann case without (with) an external potential). We will 
not use this way in the present work. 



Coming back to ( |35[) and introducing the matrix R: 

R^^l + N^l-N)- 1 (36) 

we get: 

det(if + 7) = I] W~ a 77T7 m det ( C " det (l " Q R ) ( 37 ) 

^ W a p det(l + R) 

Q = (^D-C)-\^D + C) (38) 

Let us remark that, for the free case (V(x) = 0), the matrix R is quite simple. Indeed, in 
that case: 

N (a/3)(pri) = <W<5/3r; COth y/jl a /3 + 

(39) 

R(af3)(n V ) = Sa V Spne~^ laf> ( 4 0) 

The matrix R couples any arc (a/3) to its time-reversed (f3a). Those considerations will show 
useful at the end of this paper. 

Now, let us show that, for permutation-invariant boundary conditions (and V(x) ^ 0), the 
spectral determinant can be expressed in terms of a vertex (V x V) matrix. 

4 Permutation-invariant boundary conditions 

The boundary conditions are now given by eqs.(|].|7]). 
4.1 With an external potential 

To compute the spectral determinant, we will proceed as before but, this time, we will consider, 
for each bond, two other independent solutions, Xa/3 and Xpai °f the equation (H + 7) x = 0. 
They are chosen, now, to satisfy the following conditions: 

c a Xaf3{a) + d a ^^-(a) = 1 (41) 
dx Q/ 3 

cpXam + dp^-iP) = (42) 

c a X/3a(oi) + 4v^(a) = (43) 
dx a p 

cpxM+dp^iPl) = 1 (44) 



3 For local coupling, C and D are block-diagonal matrices. N is always block-diagonal but not built with the 
same blocks as C and D. Thus, no further simplification appears in that case when evaluating det(C + DN) in 
eqM). 
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As before, we will set: 

The wronskian of Xo/3 and X/3o writes: 



W Q ^x.^-X^^=W^ (45) 
With eqs. (|4T])-(|4^) , we get the useful relations: 

c a W af 3 = ^j^-(a) ; d a W a (3 = -X/3a(a) (46) 
dx aj 3 

Let us show what happens for the Green's function G(x,y). We still assume y G [ab]. 
For x G [a/5] / [aft], we write: 

G(x, y) = B {a/3) (y)x a p(x) + B { p a) (y) X/3a (x) (47) 
where the quantities -B( Q/ g) (y) are to be determined. 

Of course, if x G [aft], an additionnal term of the form "XabXba/Wab" must appear (see eqs. 



Nevertheless, with the boundary conditions (^Q), it can be shown that, for any vertex a, 
the quantity Br a pA(y) where /3j is a nearest neighbour of a, does not depend on i. In those 
conditions, we can set: Bf a p.\(y) = B a {y) and write for the Green's function: 

i) x G [a(3] j=- [ab] 

G{x,y) = B a (y) Xa p(x) + Bp(y) X p a (x) (48) 

ii) x G [aft] 

x<y G(x,y) = £? Q (y)Xafc(x) + B b (y) Xba (x) + (49) 

VVab 

x>y G(x,y) = B a (y) Xab (x) + B b (y) Xba (x) + * 6 °^* ab ^ (50) 

The boundary conditions lead to the equation: 

MB = C (51) 
where M is a (V x V) matrix with elements: 

E XaA («) J + w « (E XU («) J ( 52 ) 

= {c a w a — t a d a )W a f3 if [a/3] is a bond (53) 
= otherwise 

B and £ are two (V x 1) vectors of components: 

B a = B a (y) (54) 

£a = ~ ($aaXab{y){CaW a ~ d a t a ) + 5 abX ba(y)(cbW b - d b t b )) (55) 



7 



Solving ( |5ll) and taking the trace of G with the relations |9| : 



d a \ xlb = diXab{a) (56) 
J [06] 



d b \ xL = d 7 Xba(b) (57) 

J [ah] 

XabXba = -d^Wba (58) 



we finally get the spectral determinant (still up to a multiplicative constant): 

det(F + 7 ) = n^ det ( M ) (59) 



Comparing the asymptotic behaviours of the right-hand sides of ( p5| ) and ( J59D when 7 —* 00, 
we establish the following equality that is valid in the presence of a potential V(;c) and for 
permutation-invariant boundary conditions: 

II det ( C + DN) = Y[ 7TT~ det(M) (60) 

[a/3] ^ [a/3] W ^ 

Recall that, for such boundary conditions, C and D are block-diagonal matrices given by eqs.(^| 
0)- 

4.2 Free case 

Let us study the case V{x) = still with permutation-invariant boundary conditions. 
With the notations 

* = (61) 



Ma 

1 + m a /ic 



(62) 



M « = (63) 



(6C, 37, p|) lead to: 



det(l - QR) = 2~ v J] (p aVa ) JJ (l - r^e" 2 ^) det M (64) 

Q [a/?] 

with the (V x V) M matrix: 

Maa = ^-^ + ^Ef 1 + W ^^ (65) 

_2e - v / 7'«/3 

M Q/3 = _ 2 jEj - if [a/3 ] is a bond (66) 

1 — VaVf3 e V' af> 

= otherwise 

For permutation-invariant boundary conditions, the matrices C, D and Q (eq.(|38|)) are 
block- diagonal. The block Q a takes the simple form: 

Qa = Va ("I + PaF a ) (67) 
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The only non- vanishing elements of the QR matrix are: 

In view of the following application, we will say that p a r] a — r\ a is the reflection factor in a 
and PaJ] a is the transmission factor. 

Expanding 

lndet(l - QR) = - ~ Tr (QR) n (69) 

n=l U 

and following the development of 0], we finally get: 



det(l - QR) = J] ( 1 - p(C)e~^ l( - c) j (70) 
c 

where the product is taken over all primitive orbits C. Recall that an orbit is said to be primitive 
if it cannot be decomposed as a repetition of any smaller orbit. 1(C) is the length of C. 

An orbit being a succession of arcs . . . (ret) (a (3) . . . with, in a, a reflection (if r = (3) or a 
transmission (if r ^ (3), the weight p(C), in eq. (|70|), will be the product of all the reflection - 
or transmission - factors along C. 

Henceforth, we will consider the situation where the spectral parameter 7 is equal to 1 and, 
in addition: 

PaVa = 1 5 Va=V ] l a /3 = I 

for all the vertices and bonds of the graph. 
With u = e~ l , (34) takes the simple form: 

H (l - (1 - rj) nR ^u m ^j = (1 - r] 2 u 2 ) B - v det ((1 - f?V)l + 7]u 2 Y - uA) (= Z" 1 ) (71) 

Cm 

m is the number of steps of the primitive orbit C m and nn(C m ) is the number of reflections 
(backtrackings) occuring along C m . 

Y is a (V x V) matrix with elements Y a n = 5 a p m a and A is the adjacency matrix (A a p = 1 
if \a(3\ is a bond, = otherwise). 



Setting rj = 1 implies n^(C m ) = in the left-hand side of (|7l|): we recover Ihara's formula [15 
16] where only primitive orbits without tails and backtrackings are kept. (Ihara |l5| established 
this formula for a regular graph; the proof for a general graph is done in [^] using a direct - 
and somewhat tedious - counting technique). 



Now, let us consider random walks with a given number of backtrackings. Eqs. 
suggest an expansion in closed random walks on the graph. Taking Z in (|7l|), we get: 

117 00 m V 

u ^=EEi: (1 - u m (72) 

m=2 p=0 a=l 

where iV^(a) is the number of m-steps closed random walks on the graph starting at a, with p 
backtrackings. 



69) and (M) 
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For the complete graph (each vertex a is linked to all the other vertices of the graph), we 
get the results: 



and also: 



7V 2 °(a) 


= 










N$(a) 


= (v- 


1)(V 


- 2) 






N%(a) 


= (v- 


1)(V 


-2)(V- 


-3) 




N° 5 (a) 


= (v- 


1)(V 


- w- 


"3)0/ 


-4) 


iV 6 °(a) 


= (v- 


1)(V 


- 2)(V 3 




+ 29V 




N%(a) = 


jvJO 


a) = Nl(a) -- 


= 




N£(a) = 




- 


-2)(y 


-3) 




N*(a) = 


6(V 


-1)(V 


-2)(V 


-3) 2 



(73) 



(74) 

In [p!^ 1, the same problem is studied with probabilistic methods but for open random walks. 
Closed walks are therefore obtained by identifying the starting and ending points but nothing 
is said about an eventual backtracking occuring at that point. So, the results of |l7j (let us call 
them jV^(a)) will, in general, differ from ours. For instance, we checked for the complete graph, 
the relationship: 

AC (a) = iV° (a) + -N l m (a) (75) 
m 

(This comes from the complete symmetry of this graph). 



5 Conclusion 

We have computed the spectral determinant for a Schrbdinger operator on a graph with quite 
general boundary conditions. The result is expressed in terms of an arc matrix. When the 
conditions are permutation- invariant, another expression can be derived in terms of a vertex 
matrix. Comparison of both expressions allowed us to study reflection properties of random 
walks on any graph. 

The expansion ( |70"|) of the spectral determinant in periodic orbits is the basis for obtaining 
a trace formula (see, for instance, § where this is done in great details for Neumann boundary 
conditions). Unfortunately, in the general case, the reflection and transmission factors are 7- 
dependent and technical difficulties prevent from getting a trace formula in an appealing form. 
So, this problem is still an open one. 



I acknowledge Pr A Comtet and Dr C Texier for stimulating discussions. 
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